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Abstract 

We study the superfluid-insulator transition of a particle-hole symmetric sys- 
tem of long-range interacting bosons in a time-dependent random potential 
in two dimensions, using the momentum-shell renormalization-group method. 
We find a new stable fixed point with non-zero values of the parameters 
representing the short- and long-range interactions and disorder when the in- 
teraction is asymptotically logarithmic. This is contrasted to the non-random 
case with a logarithmic interaction, where the transition is argued to be first- 
order, and to the 1/r Coulomb interaction case, where either a first-order 
transition or an XY-like transition is possible depending on the parameters. 
We propose that our model may be relevant in studying the vortex liquid- 
vortex glass transition of interacting vortex lines in point-disordered type-II 
superconductors. 
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In recent years, there has been great interest in quantum critical phenomena occurring 
in systems of interacting bosons ]!]-§]. This has largely been motivated by several beautiful 
experiments on the superconductor to insulator transition in thin metallic films at very low 
temperature [P-|T4|j. In many models, the superconductor (or charged superfluid) to insu- 
lator (or Bose glass) transition in disordered films is considered as a localization transition 
of bosonic Cooper pairs in an external random potential In non-random systems, 

the superfluid to (Mott) insulator transition can be triggered by changing the density of 
bosons through critical values commensurate with the underlying periodic lattice |l|J^j . Since 
Cooper pairs are charged, one needs to consider the influence of long-range interactions on 
the nature of the critical phenomena. In the absence of disorder, the superconductor-Mott 
insulator transition of charged bosons in two dimensions interacting by the 1/r Coulomb 
potential has been studied by Fisher and Grinstein via renormalization-group methods [|TJ. 
They showed that when the (bare) Coulomb interaction is sufficiently weak, the transition 
belongs to the three-dimensional XY-model universality class, whereas when the Coulomb 
interaction is sufficiently strong, it is first-order. In this paper, we generalize Fisher and 
Grinstein's model to the cases with the 1/r or logarithmic interaction in a time- dependent 
random potential. Time dependence of the potential in classical static critical phenomena is 
unimportant |§. In quantum critical phenomena, however, it can lead to new universality 
classes as will be shown in the present work. As well as being relevant in studying novel 
quantum critical phenomena, our model may be also a useful representation of a system 
of long-range interacting vortex lines in point- disordered type-II superconductors in an ex- 
ternal magnetic field. According to the mapping originally due to Nelson, vortex lines in 
three spatial dimensions are considered as the world lines of bosons residing in two spatial 



dimensions perpendicular to the magnetic field [15HT7|. The dimension parallel to the mag- 
netic field is mapped to the imaginary time dimension of a bosonic system. Therefore the 
static random potential provided by isotropically distributed impurities in superconductors 
can be regarded as a time-dependent potential that is random in both space and time. The 
quantum superconductor-insulator transition of interacting bosons is analogous to the clas- 



sical vortex liquid-vortex glass transition of interacting vortex lines. We caution that in our 
model, the long-range interaction between vortex lines is confined to the plane perpendic- 
ular to the magnetic field. Applying the standard momentum- shell renormalization-group 
method to the logarithmically interacting case, we discover a new stable fixed point in two 
dimensions with non-zero values in the short- and long-range interactions and disorder when 
there exists a time-dependent random potential, whereas we find a first-order transition in 
the non-random case. Since the interaction between vortex lines is logarithmic in a wide 
range of length scales when the magnetic field is close to the upper critical field [fL8|l , we 
propose a possibility that our new random fixed point may describe the real vortex liquid- 
vortex glass transition in point-disordered type-II superconductors. On the other hand, 
when the interaction is the 1/r Coulomb interaction, we find a behavior analogous to Fisher 
and Grinstein's results in a non-random system. 

We consider the (d + l)-dimensional classical action for the m-component classical com- 
plex Bose field ijji(x,r) coupled to the scalar gauge field A(x, r) mediating the long-range 
interaction between charge e bosons in d spatial and one temporal dimensions: 
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+ \j:(ck°+gk 2 )\A(k,u;)\\ (1) 

where A(k, to) is the Fourier transform of A(x, r). a, b, u, c and g are positive constants and 
r(x, t) is a Gaussian random function satisfying 

r(x, t) = r + Sr(x, r), (8r(x, r)) = 0, 

(5r(x, T)5r{x', r')) = 2v8(x - x' )6{t - r), (2) 

where (• • •) denotes the statistical averaging and v (> 0) measures the strength of random- 
ness. The action S applies to the cases where there is a particle-hole symmetry and the 
interaction potential between bosons is proportional to \jr d ~ u in the long distance limit, r 



being the distance between bosons. The 1/r Coulomb interaction corresponds to a = d— 1, 
while the asymptotically logarithmic interaction corresponds to a = d. In the absence of the 
particle-hole symmetry, one has to replace the term proportional to a in (|1|) by a term of the 
form ipi*(d T — ieA)ipi 0,0]. In general, disordered systems do not possess the particle-hole 
symmetry. We note, however, that it has been argued recently, in a rather compelling man- 
ner, that the particle-hole asymmetric term is irrelevant in the renormalization group sense 
and the exact (statistical) particle-hole symmetry is restored after averaging over disorder 
at the critical points of various kinds of disorder- driven superfluid transitions, including the 
vortex glass transition Hl9|-|2^| . We accept these arguments and assume that the particle- hole 



symmetric action, Eq. (|l|), describes the critical fixed point of the vortex glass transition. 
When r is non-random and g is zero, the superconductor to Mott insulator transition de- 
scribed by this model was studied by Fisher and Grinstein |]]. In the present paper, we 
generalize their work to the cases where r is a Gaussian random function of both space and 
time. 

We use the standard replica trick to reduce S to a non-random effective action: 

p m 
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«,7=1 i,j=l 

+ \tE(^ + 9k 2 )\A a (k,uj)\\ (3) 
Q=1 k,u 

where a and 7 are replica indices and the number of replica p is sent to zero at the end of 
the calculation. In order to apply the momentum shell renormalization group method, we 
perform a Fourier transform of Eq. (|3|). The Fourier-transformed action has five types of 
interaction vertices. We perform a perturbation expansion in terms of these vertices to one 
loop order and integrate out the components with wave number k in a shell k\/e l < k < k\, 
where k\ is the momentum space cutoff and e l is the rescaling parameter. We recover the 



original form of the action by renormalizing r , a, b, g, u and v using the results of the 
momentum space integration and rescaling frequency, momentum and the fields with the 
rules 



where z is the dynamical exponent and and (a are rescaling factors of the fields. 

We perform our calculations in an e (= 3 — d) expansion about the upper critical dimen- 
sion d = 3. Firstly, we consider the competition of the k a and k 2 terms in the gauge field 
propagator. A naive power counting suggests that the k 2 term is irrelevant as long as the 
k a term with < a < 2 is present. This conclusion is not always true, however, since the k 2 
term is renormalized in a perturbation theory, whereas the nonanalytic k a term is not, and 
in some cases, this renormalization makes the k 2 term more relevant than the k u term. To 
confirm this, we fix a = & = g = 1 and perform a perturbation calculation assuming that c 
is irrelevant. To one loop order, we obtain the exponents 



k^k/e l , 
A a (k,u) - 



uj — > uj/e , 
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(4) 



z 




2Q = l-d+- 



(5) 



and the recursion relations for c: 




(6) 



and for r , u and v: 



where u = UKd/A, v = VK4/4, w = wkj = e 2 /^ and Kd = 2/ (4:7i) d ^ 2 T(d/2). The momentum 
cutoff was chosen to be equal to f. The recursion relation for the parameter w = e 2 is 
obtained from the condition that the action is invariant under the gauge transformation 

= ^ e te \ A = A + ^ (8) 
dr 

for any function A(r) to all orders. We find 

— = (2Ca + 2z)w = ew- J w 2 , (9) 

which possesses a stable fixed point w* = 12e/ (m + 4). Substituting this value into Eq. @, 
we observe that the parameter c is indeed irrelevant at this fixed point if 2 > o > o c = 
2 — mw* : /12 = 2 — me/ (m + 4). In the cases with < a < a c , where g is irrelevant, we fix 
a = b = c = 1 and g = and perform a one-loop calculation. It turns out that the recursion 
relations fl7|) remain the same, whereas Eq. (|9|) is replaced by 

_ = ( f - £> - T , (10) 

where we have introduced another small parameter e CT = 2 — a. 

We consider the 1/r Coulomb interaction case first. According to the results obtained 
above, the k 2 term in the gauge field propagator is irrelevant to the first order in e in the 
physically interesting case with e = 1 (that is, d = 2) and m = 1. Then e — t a = 1 — d+cr = 
and the marginally irrelevant parameter w satisfies 

f-f- w 

Integrating this equation and the recursion relations for u and v numerically for various 
initial values of u, v and w, we find the renormalization group flow diagram qualitatively 
similar to that discovered by Fisher and Grinstein. That is, there is a plane in the parameter 
space separating the domain of attraction to the (4 — e)-dimensional random XY-model fixed 
point with positive u* and v* and zero w* [53| and the domain where all points flow to the 



region with large negative u, large positive v and zero w. When the physical value of w is 
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sufficiently small, the flows are attracted to the XY fixed point. When the initial w value is 
large, u grows to a large negative value, v grows to a large positive value and the effective 
action (|3[) becomes thermodynamically unstable. This behavior is usually interpreted as a 
signature of a fluctuation-driven first-order phase transition. We illustrate this behavior by 
showing the flow diagram in the non-random case with v = in Fig. 1. 

In case of the logarithmic interaction in two dimensions, we have only the k 2 term in the 
gauge field propagator. In this case, the effective charge parameter w has a stable fixed point 
at w* = 12e/ (m+4). Let us consider the non-random case with v = first. Then it is easy to 
show that u possesses a stable fixed point at u* = (m + 10 + \Jm 2 — 52m — 188)e/4(m + 2) 2 
only when m > m c = 26 + 12\/6 ~ 55.394. When m is smaller than m c , all points flow to 
large negative u values as shown in Fig. 2. We cautiously interpret this as a signature of a 



first-order phase transition p4| . In the random case, we substitute the fixed point value of 
w into the recursion relations for u and v and find a stable fixed point with positive u*, v* 
and w* when m — 1. The fixed point values are 

ll + v/409 ^ v/409 ^ 12 /io . 

u = e, v* = e, w* = — e. (12) 

40 40 5 y J 



This point has three irrelevant eigenvalues — e and —(33 ± -\/2183i)e/20. The complex 
eigenvalues suggest that the renormalization-group flows spiral into the stable fixed point 
and give oscillatory corrections to scaling. In Fig. 3, we show the flow diagram in the (u, v) 
plane with w fixed at w*. There is a separatrix dividing the domain of attraction to the 
random fixed point and the unstable region. The critical exponents associated with the new 
random fixed point are 

2 ,4 1 n ( 9 v / 409\ 

ri = -e, z — 1 e, — = 2 — 1 e, 

' 5 ' 5 ' v V 10 20 ) ' 



i = l = 2+a-^U (13) 
v r zv V10 20 J ' y ' 

where v and v T are the correlation length and time exponents and r] is defined by 2^+d+z = 
2-rj. 
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In conclusion, we performed a renormalization-group analysis of the superfluid-insulator 
transition in a particle-hole symmetric model of long-range interacting bosons with or with- 
out a time-dependent random potential. We find a new stable fixed point with non-zero 
disorder and long-range interaction when the interaction is logarithmic. Without a random 
potential, the transition in logarithmically interacting systems is argued to be first-order. We 
propose a possibility that our model may be relevant in interpreting the vortex liquid- vortex 
glass transition of interacting vortex lines in type-II superconductors. 
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FIGURES 




FIG. 1. Flow diagram for the non-random Bose system interacting by the 1/r Coulomb poten- 
tial when m = 1. S labels the separatrix dividing the domain of attraction to the XY fixed point 
and the unstable region. 
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FIG. 2. Flow diagram for the logarithmically- interacting non-random Bose system when m = 1. 
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FIG. 3. Flow diagram for the logarithmically-interacting Bose system in a time-dependent 
random potential when m = 1 and w = w* = 12e/5. S labels the separatrix dividing the domain 
of attraction to the random fixed point with positive u, v and w values and the unstable region. 
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